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ABSTRACT 

We have applied the SRIM computer code to study the sputtering of some likely 
astrophysical grain materials, and we have shown that selective embedding of metallic 
projectiles offers a partial explanation of gas-phase depletions. We show that supernova 
Shockwaves sweep a significantly larger mass of interstellar gas per unit time than the 
Shockwaves generated by outflows in star-forming regions. We apply our sputtering 
model to the bombardment levels expected in a supernova shock, and show that net 
embedding may dominate over net sputtering, leading to grain growth under some 
circumstances, particularly when the bombarding gas is enriched with metals from the 
supernova progenitor star. A combination of short cooling times and net embedding 
mean that it is possible for a type II supernova to generate more dust that it destroys, 
and we conclude that, in general, the sputtering process often leads to a compositional 
change in the grain material rather than simply to grain erosion. 

Key words: ISM: dust — ISM: evolution — atomic processes 



1 INTRODUCTION 



The life-cycle of astrophysical dust is not well understood. 
Whilst certain constraints can be set on the class of stars 
which are sources of dust (Whittet (1992), Jones (1997)) and 
additional information can be gleaned from isotopic compo- 
sitions of grains (Zinner 1998) we know surprisingly little 
about the fate of an average grain from the point where it 
passes from the stellar wind of its host star, or supernova, 
into the wider interstellar medium (ISM). 

Of particular interest for models of the life-cycle of dust 
(Edmunds 2001) are a set of competing processes which con- 
trol the size and mass of the grain, some tending to increase 
the average grain mass and others degrading or even de- 
stroying the grain. The effect of many of these processes on 
grains were considered by Barlow & Silk (1977) and Barlow 
(1978a,b). One way to divide up these processes is to classify 
them cis single particle processes, where a single gas atom or 
molecule collides with the grain, and grain-grain processes 
where there is a collision between two lumps of solid ma- 
terial. In both cases the result of the collision, in terms of 
accretion or erosion of the grain, depends on the details of 
the collision, such as the energy, angle, and composition of 
the coUision partners. Such processes have been included in 
sophisticated models of chemical evolution, which include 
the development of Galactic dust properties, for example 
Dwek (1998) and references therein. 

Single particle collisions with a grain may result in ac- 
cretion of gaseous material onto the grain, or in sputtering 



when the incoming particle has sufficient energy to eject a 
grain particle into the gas. In cold molecular clouds, we can 
deduce that the single particle processes lead to accretion of 
gas molecules onto the grains, leading to an increase in the 
average grain size, and mass (de Muizon et al. 1998). The 
key observations in this respect are signatures of ices on dust 
in the spectra of dark clouds (Whittet et al. 1996). More- 
over, we can see that as the molecular cloud regions become 
colder, more and more volatile ices are deposited on the 
grains, with signatures of non-polar ices, such as non-polar 
CO, appearing instead of the water-ice and other polar ices 
common in higher temperature regions (Gibb et al. 2000). 

Far less is known about the consequences of single par- 
ticle collisions in the warmer and more diffuse ISM which 
surrounds dark clouds. Observations are sufficient to show 
that in these regions, grains have lost their ice mantles, ex- 
posing bare core material. There is evidence (Whittet 1992) 
that this is often of silicate composition, but may also be car- 
bonaceous. It is also quite possible that refractory organic 
material acts as a kind of inner mantle (Greenberg & Li 
1996) protecting a true inner core of silicate. The absence of 
ice mantles is often taken to show that grain erosion by sput- 
tering must be an efficient process, particularly in the warm 
ISM. Although the thermal energy of gas atoms in this phase 
is typically well below the surface energies of likely grain 
materials, making thermal erosion impossible, thermal and 
non-thermal sputtering in supernova shocks is usually con- 
sidered to be adequate to supply sufficient gas or plasma par- 
ticles with the necessary collision energies (McKee (1989), 



© 0000 RAS 



2 M. D. Gray & M. G. Edmunds 



McKee et al. (1987)). Indeed, supernova shocks are taken by 
many authors to be the main sites of grain destruction in the 
ISM, for example, Itoh (1985), McKee (1989). The extent 
of destructive processes is, however, rather controversial, 
since there is evidence from SN1987A (Kozasa, Hascgawa & 
Nomoto 1991), (Arendt, Dwek & Musely 1999) and Cas A 
(Dunne et al. 2003) that supernovae can also produce new 
dust - although its subsequent fate is unknown. Theoretical 
condensation models (Todini & Ferrara 2001) and general 
consideration of grain abundances (Edmunds 2001) would 
argue for supernovae being a significant source of interstel- 
lar grains, implying that destruction in these environments 
cannot dominate. 

Obviously the main single particle collision partners for 
grains will be hydrogen and helium, and in supernova shocks 
these will probably be in the form of ions. Heavier ions, 
which carry more momentum, are more effective at sput- 
tering than their abundances would indicate. However, an 
extensive study (Field et al. 1997) indicates that their to- 
tal effects are probably not more than those of the hydro- 
gen and helium, at least for typical Galactic abundances of 
metals and an energy spectrum of projectiles which extends 
well above the sputtering threshold (sec Section 3.1). More 
recent work (May et al. 2000) has studied the effects of sput- 
tering of typical silicate core materials in considerable detail, 
though restricting the bombarding particles to helium and 
heavier species. In this work, we present a detailed account 
of the sputtering process, via the Monte-Carlo ion impact 
computer code SRIM (Ziegler, Biersack & Littmark 1985), as 
in (May et al. 2000), but consider phenomena not covered 
by that work: timoscales for the release of grain material, 
sputtering by hydrogen, sputtering of ice-mantles, and the 
study of sputtering as a process of grain modification, rather 
then the simple errosive process which is normally assumed. 



2 THE SRIM PACKAGE AND PARAMETERS 

The SRIM computer package was originally written to solve 
problems in nuclear physics, where the penetration of a tar- 
get by ionizing particles is to be investigated. The program 
takes a Monte-Carlo approach, in which the tracks of in- 
dividual particles through the target are followed, and the 
results from many such independent tracks are accumulated 
to derive useful statistics for the target/projectile combina- 
tion. It should be emphasized that the individual tracks are 
completely independent, and assume a virgin target: that is 
projectiles which are stopped within the target are assumed 
not to change its composition. Part of this work (see Sec- 
tion 3.3 and 3.4) lifts this assumption by following a set of 
calculations, each of which has a target composition modi- 
fied according to the results of the previous computation. 

The basic physics contained in the SRIM package relates 
to various nuclear stopping processes, by which the input 
projectile loses its original energy to particles in the target. 
Secondary processes, or cascades, are dealt with fully, that 
is any target particle which recoils with a significant fraction 
of the energy of the projectile is then treated as a further 
projectile which then also has to lose its energy to further 
target atoms, or escape from the target. There are various 
parameters involved in running SRIM: the target composition 
and thickness, projectile type and energy, projectile angle. 



and various operational flags which control the form of the 
output, and the detail in which the projectile tracks are re- 
ported. There are also three energies related to the lattice of 
the target: the surface energy, Es, the displacement energy, 
Ed, and the lattice binding energy, Eb- The first of these is 
the energy needed to eject a surface atom from the target. 
The second is the energy required to displace an atom from 
its lattice site, measured from the bottom of the potential 
well in which the target atom can move, to the top of the 
potential barrier that nmst be overcome to leave the site. 
Once free of its lattice site, a target atom may be able to 
escape the target as a sputtered atom, but if it does escape 
to infinity, it must have gained an energy at least equal to 
the lattice binding energy, relative to when it was bound in 
its lattice site. 

Most of the work reported here involves the bombard- 
ment, by protons and heavier ions, of targets which approx- 
imate reasonably to interstellar dust surfaces. We note that 
SRIM has already been used to study the effects of sputtering 
by heavy ions in the CNO group (Field et al. 1997), and for 
more limited results with hydrogen and helium (Tielens et 
al. 1994). The versions of the SRiM code used in this work 
were SRlM-2000.10, for the comparison with Tielens et al. 
(1994), and and SRIM-2000.40, for the work on olivine sub- 
strates. 



3 SPUTTERING RESULTS 

The results of this work fall into the following categories: 

First, we verify the SRIM method, by comparing our sput- 
tering yields, as a function of projectile energy and angle, 
with the work of other authors, both for SRIM calculations. 
Field et al. (1997), May et al. (2000), and for other methods 
(Tielens et al. 1994). Secondly, we consider the effects of se- 
lective sputtering by introducing small amounts of metallic 
impurities into the grain material; we discuss whether or not 
this could help to explain the observed gas-pliasc dcflicien- 
cies in certain elements. Thirdly, we discuss the overall effect 
of the bombardment of grain materials by common ions, in- 
cluding the modification of the grain material by projectiles 
which are stopped inside the grain. 

3.1 Comparison with Previous Work 

In Figure 1, we plot the sputtering yield, Y{E), defined as 
the number of target atoms ejected per projectile, expressed 
as a percentage , as a function of projectile energy, for three 
different target materials. The projectiles are taken to strike 
the target at an angle, 9, of zero with respect to the surface 
normal. The projectiles are hydrogen nuclei in all cases, and 
the energy range is from 2 eV to 2 MeV. The target materi- 
als have been chosen to represent three likely populations of 
astrophysical grains: graphite and silicate to represent grain 
cores, and water-ice to represent grains which have accumu- 
lated cold mantles. In Figure 2, for a more restricted range 
of projectile energy, we show the angle- averaged sputtering 
yield, assuming an isotropic distribution of projectile trajec- 
tories. This average sputtering yield is given by 

Y{E) = / Y{E,9)sm9de (1) 
Jo 
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Figure 1. Percentage sputtering yields for liydrogen nuclei, nor- 
mally incident upon surfaces of graphite, silicate, and water-ice, 
as a function of the projectile kinetic energy. 
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Figure 2. Angle-averaged percentage sputtering yields for hydro- 
gen nuclei, incident upon surfaces of graphite, silicate, and water- 
ice, as a function of the projectile kinetic energy. An isotropic 
distribution of projectile trajectories has been assumed. 



The main features of Figure 1 are the following: Sput- 
tering is an essentially inefficient process; we see that even 
at the optimum energy for the most susceptible material 
(ice) the yield is less then 10%. We also see that there is 
a band of interesting energies. For the target used here, of 
thickness 100 nm, we find little sputtering outside the range 
30 eV to 100 keV. The lower figure depends only on the sur- 
face energy, displacement energy and lattice binding energy 
of the material. On the other hand, the upper limit depends 
on the grain size, so that we expect significant sputtering 
damage to persist to higher energies for large grains, whilst 
very small grains will become immune to sputtering at lower 
energies than in Figure 1. This view of the sputtering pro- 
cess is reinforced by Figure 2: even with the bias to large 
angles introduced by the solid-angle integral in eq.(l), no 
mean sputtering yield exceeds 20%. 



Field ot al. (1997) and May ct al. (2000) do not con- 
sider sputtering by hydrogen projectiles, so we compare our 
results in Figure 1 and Figure 2 to the work of Tielens 
et al. (1994), which also contains some experimental data. 
For carbon (graphite) the top left-hand graph of Figure 10 
of Tielens et al. shows hydrogen sputtering from a carbon 
(graphite) surface. The peak yield is about 1% at 350 eV, 
and the threshold for sputtering is axound 40 eV. In our work 
(Figure 1) we find a peak yield of about 1% for impact en- 
ergies near 1 kcV, and a sputtering threshold energy near 
50 eV. The hydrogen sputtering in Figure 1 of Jurac, John- 
son & Donn (1998) has a similar peak yield, but is shifted to 
slightly lower energy than the equivalent graph in Tielens et 
al. The silicate graph in Figure 10 of Tielens et al. docs not 
have a plot for hydrogen, but comparing it with the graph 
for SiC, and assuming the same order of magnitude differ- 
ence between H and He sputtering, our figure of about 3% 
for the peak sputtering yield from hydrogen impacts seems 
reasonable. For ice, we find a larger discrepancy between our 
results and those of Tielens et al.: we calculate a peak sput- 
tering yield for hydrogen ions, normally incident on ice (see 
Figure 1), of 6.5%, but Tielens et al. find 30%. It is not clear 
where this difference arises, but we note that Tielens et al. 
have scaled their theoretical curves via the free parameter, 
K, in order to fit experimental data, of which only a small 
amount is available. One likely explanation is that neutral- 
ization of a projectile ion, when striking an ice surface, can 
lead to electronic sputtering, a process not included in SRIM 
(Jurac, Johnson & Donn 1998). Another is that that there 
are many allotropes of ice, each with its own lattice bind- 
ing energy. However, as we do not consider ice further in 
this work, but concentrate on silicate cores, the discrepancy 
noted above is not particularly important here. 

For comparison with May et al. (2000) we choose an 
olivine target, and look at the sputtering yield of silicon 
when the olivine, MgFeSi04, is bombarded by oxygen and 
iron projectiles (see Figure 2 of May et al). When we reverted 
our SRIM code to include the restrictions imposed by their 
earlier version, the results were the same to within the ac- 
curacy with which we could read their Figure 2 (~ 5%). For 
example, at 100 eV with oxygen projectiles, we obtained a 
yield of Si of 0.018; the marker on Figure 2a of May et al. 
at 100 eV is clearly just below 0.02. For the same graph, 
at 60 eV, we find a yield of Si of 0.00204, whilst the graph 
value at this energy is almost exactly 0.002. This excellent 
agreement was found in spite of a small difference in the 
density used for the olivine (we used 3.81 gcm~^ (The CRC 
Handbook 1996) instead of the 3.84 gcm~^ used by May et 
al.), although the yield is expected to be quite insensitive 
to the density: Field et al. (1997) state that a 5% change 
in sputtering yield requires a change in density as large as 
20%. Note that we have used the fractional definition of the 
yield here, eis in May et al., rather than the percentage defi- 
nition used elsewhere in this work. The principal restriction 
of the older version of SRIM used by May et al. is that it 
allowed only a single value for each of the energy parame- 
ters, Es, Ed and Eb, defined in Section 2. Therefore, they 
used a weighted average over all the constituents of each 
target. In our work, we used the more recent SRlM-2000.40 
to study the olivine target, and this allows separate values 
of Eg, Ed and Eb to be chosen for all species in the target. 
We used the values supplied by SRiM-2000.40 which appear 
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Table 1. The Energy Parameters Eg, E^i and Eg for Olivine 



Element 


Es 


Ed 


Eg 




eV 


eV 


eV 


Mg 


1.54 


25.0 


3.0 


Fe 


4.34 


25.0 


3.0 


Si 


4.7 


15.0 


2.0 


O 


2.0 


28.0 


3.0 



Sputtering iiomGiaphite 



in Table 1 above. None of the individual species energies 
supplied by SRIM are as high as the blanket averages used 
in May ct al. It is therefore not surprising that, in our work 
on olivine, we find significant sputtering of Si by O, Fe, and 
other projectiles, at 30 eV. When using the species-specific 
energy parameters, we obtain a fractional yield of silicon, at 
100 eV, of 0.055 for the oxygen projectile and 0.034 for the 
iron projectile. Below this energy, large dcscrcpancies appear 
between our data and May et al. (2000). Fortunately, we are 
considering a higher energy regime than May et al., so the 
precise energy of the sputtering threshold is less important. 

In Figure 3, Figure 4, and Figure 5 we show more detail 
in the sputtering yield, including, in these figures, the de- 
pendence on the angle of incidence of the projectile. Figure 3 
shows the results for a carbonaceous grain (graphite). Fig- 
ure 4 shows analogous results for a silicate material (quartz), 
and data for water-ice are depicted in Figure 5. All the sur- 
faces exhibit a qualitatively similar behaviour: higher angles 
of incidence increase the sputtering yield for a given energy, 
until an optimum angle, usually in the range 75-88 degrees. 
At higher angles, the number of backscattcred projectiles bo- 
comes so large that sputtering efficiency is lost. Only near 
the optimum energies for the most susceptible material (ice) 
does the sputtering yield at the optimum angle exceed 30%. 
The form of the angular dependence in Figure 3 to Figure 5 
can be compared to Figure 2 of Jurac et al. (1998), noting 
that the latter results are for a helium projectile. 

A conclusion that we share with Tielens et al. (1994) is 
that thermal sputtering is rarely of importance in the qui- 
escent ISM. Only the hot, or coronal gas, phase of the ISM 
has sufficient temperature to place ions above the sputtering 
threshold of typical grain materials; this hot gas is, however, 
too diffuse to erode grains within a reasonable time. Sputter- 
ing therefore proceeds rapidly only in shocks, in which some 
of the particles at least are subject to non-thermal accelera- 
tion mechanisms and/or thermal sputtering at much higher 
densities than exist in the coronal gas. The two likely can- 
didates for shock-sputtering of grains are supernova shocks, 
for example Tielens (1994), Jenkins (1989), McKee (1989), 
and the C-type shocks fronting molecular outflows in star- 
forming regions, as considered by Field et al. (1997) and May 
et al. (2000). Which of these shock processes is dominant? 

Thermal sputtering is out of the question in cold molec- 
ular clouds. Therefore, Field et al. (1997) and May et al. 
(2000) introduce non-thermal sputtering via a model of 
a C-type (continuous) MHD shock, of modest speed (20- 
45kms~^), moving through a molecular medium. The post- 
shock gas does not reach temperatures sufficient to ionize 
most components of the gas, so neutrals remain abundant, 
and sputtering results from a large drift velocity between 
the neutral and ionized fluids. The grains are talcen to be 
part of the latter fluid; see the discussion in Weingartner 



Y(%) 




Angle (deg) 



Figure 3. Percentage sputtering yields for hydrogen nuclei inci- 
dent upon graphite, as a function of the projectile kinetic energy 
and angle of incidence. An angle of 90 degrees is parallel to the 
target surface. 



& Draine (2001) for the range of possible charge-states that 
are accessible to astrophysical dust grains: both positive and 
negative charging is possible. We now calculate the likely 
sweeping rate of the ISM by shocks of this type, assuming 
that they are produced by bipolar outflows, associated with 
star- formation. Bontemps et al. (1996) derive an expression 
for the 'momentum flux' in an outflow, based upon obser- 
vations. This flux is related to the bolometric luminosity of 
the source, which is, in turn, assumed to be produced by 
accretion during the outflow phase. Combining expressions 
for the momentum flux and accretion power from Bontemps 
et al. yields 



F{M) 



= 10-'-® log ( 

Y 3RqLq 



(2) 



where the 'momentum flux', F{M), is measured in units of 

solar masses per year times an outflow line- width in kms~^, 
and Mace is the accretion rate. Assuming a timescale for 
the outflow phase of to = rs x 10^ yr (Saraceno et al. 1996), 
where tb is the outflow lifetime in units of 10^ yr, and a veloc- 
ity width of V50 X 50kms~^, the rate at which ISM material 
is shocked by the outflow is given by 10'^Tr,F{M)/vrM, where 
V50 is the velocity width in units of 50kms^^, or 



M(M) = 



10" 



50t;B0 V 9-46 x W^'^RqtsLq 



Mqvt-^ (3) 



where we have assumed that the star accretes its own mass 
during the time To, so that Mace = M/tq. We note that 
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Figure 4. Percentage sputtering yields for hydrogen nuclei inci- 
dent upon crystalline silicate (quartz), as a function of the projec- 
tile kinetic energy and angle of incidence. An angle of 90 degrees 
is parallel to the target surface. 



Figure 5. Percentage sputtering yields for hydrogen nuclei inci- 
dent upon water-ice, as a function of the projectile kinetic energy 
and angle of incidence. An angle of 90 degrees is parallel to the 
target surface. 



vo is essentially twice the shock velocity of the outflow. As- 
suming a modified power-law initial mass function (IMF), 
with an index of 1.35 for stars below 0.5 M0, and the classi- 
cal Salpeter index of 2.35 for stars of greater mass (Kitchin 
1987), we find that, for stellar mass limits of 0.08 Mq, and 
100 M0, the mean stellar mass is 0.5253 Mq and the mean 
sweeping rate, found from 



M 



2 X lO"'' *' 



W50 



9.46 X IO^^tsRqLq 



■0.5Mq 

0.2367Mq^^ ' '""■'^^ 



M'^-'^'dM 



+ 0.1184Af, 



© 



O.OSMq 
IOOMq 



M-°-^^dM 



M@yr-' (4) 



is 5.78 X lO-V(i;5oT-5°-^)M0yr-i. We note that if the IMF 
were extended to brown-dwarf masses, the mean stellar 
mass, and mean sweep-rate would be lower than calculated 
here. Finally, the amount of ISM swept by the outflow of an 
average star, during its formation, is found by multiplying 
the sweeping rate by tq, giving 



M„„t = 0.578r5''-Vt'50 Mq 



(5) 



where we note that the dependence on the outflow lifetime 
is very weak. 

We compare the swept mass of ISM from eq.(5) to a 
figure for supernovae, derived by McKee (1989), of 6800 M© 



before the shock decelerates to lOOkms"^. Using the same 
IMF as used above for the outflows, we flnd that the frac- 
tion of stars formed as supernova progenitors (M > 8Mq) is 
0.0051. Assuming that the modern supernova rate is equal 
to the formation rate of such stars, then the relative effi- 
ciency of shocking interstellar gas by supernovae and out- 
flows is in the ratio 0.0051 x 6800 x V50 /{0.578t° '^), which is 
equal to 600v5o/t5 '^. Supernovae therefore shock about 600 
times more material than star-forming outflows, so whilst 
the outflows can be locally very important, most sputtering 
occurs in supernova shocks unless the sputtering efficiency 
in outflows is vastly more efficient than in supernovae; this 
is unlikely to be the case. 

As a check on the above calculation, we compare ob- 
servational estimates of the current supernova and star- 
formation rates, which are independent of any assumed IMF. 
McKee estimates the effective supernova rate (for super- 
novae which interact strongly with the ISM, which effec- 
tively excludes type la supernovae, and allows for correla- 
tions due to multiple supernovae in clusters) to be about 
1 event per 120 yr. We compare this with estimates of the 
Galactic star- formation rate by Boissier & Prantzos (1999), 
who estimate a rate of 3 x 10~® Mq yr~^ pc~^ for the solar 
neighbourhood. From Figure 2 of Boissier & Prantzos, we 
can see that the star-formation rate does not reach ten times 
this value at any radius in the Galactic disc. Adopting a 
Galactic radius of 15kpc, and using the mean stellar mass, 
calculated earlier, of 0.5253 Mq, we find a star-formation 
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rate of 4.0 yr^^. The ratio of the effective supernova rate to 
the total star-formation rate is therefore 0.002, compared 
with the 0.0051 calculated on the basis of our IMF. These 
figures are consistent if we adopt the lower figure, because 
the value of 0.0051 is not corrected for the effectiveness of 
the supernovae. With the lower figure, the supernovae are 
still about 240 times as effective at shocking interstellar gas 
then star-forming outflows. 

In the light of the analysis above, the model we adopt 
in the present work is very different from May ct al. (2000). 
We consider a J-type supernova blast-wave moving at over 
lOOkms"^, which provides both thermal and non-thermal 
sputtering projectiles. The post-shock gas is ionized by the 
supernova shock, so the projectiles and grains are both part 
of the charged fluid, and the iou-ncutral drift velocity mech- 
anism is not applicable. In the present work, we ignore beta- 
tron acceleration of the grains and projectiles, and also the 
fact that a grain will initially have a large initial velocity 
relative to the gas until its motion is thermalized. We divide 
the post-shock gas into a thermal fraction, which has an en- 
ergy distribution function given by the Maxwell-Boltzmann 
form. 



p{E) 



7rV2 (fcT2)3/2 



-E/(kT2) 



(6) 



rather than a skewed Maxwellian (used, for example, in Bar- 
low (1978a)), and a small non-thermal fraction, which enters 
a flrst-order Fermi acceleration process above some injection 
energy, Ei^j . We use the results of a diffusive shock model 
(Berezhko & Ellison 1999) to obtain our non-thermal energy 
distribution functions and the post-shock conditions. The 
crucial difference between these models and simpler test- 
particle calculations is that penetration of the pre-shock gas 
by energetic particles leads to energy loss from the shock: 
there is no longer a simple limiting compression ratio for the 
shock, and it can become very large. However, for acceler- 
ated particles which have not reached relativistic energies, 
the power-law spectrum of the non-thermal distribution is 
not different to the test-particle case (Berezhko & Ellison 
1999). The power-law energy spectrum of ions of species j, 
derived from Berezhko & Ellison, is 



^3/4^3/4^3/4 



23/44r,,, 



\m ) 



3/4 



E 



-7/4 



(7) 



where rj is the injection fraction of particles into the non- 
thermal mechanism, A is the ratio of the injection momen- 
tum to the thermal momentum mHCs2, "ig is the ratio of 
specific heats in the thermal gas, fh is the mean particle 
mass, and rrij is the mass of species j. The temperature in 
the post-shock gas is T2, and the adiabatic sound speed in 
the post-shock gas is Cs2 . The power-law index results from 
taking the compression- ratio in the sub-shock, rsut, to be 3 
(see Berezhko & Ellison (1999)). A crucial quantity is the 
post-shock temperature, T2. We adopt the most physically 
realistic case which allows for loss of energy from the shock, 
allowing the overall compression ratio to rise above its clas- 
sical limit of 4, but we allow the compression ratio to be 
limited by the dissipation of Alfven waves, so it cannot rise 
without limit as the Mach number tends to infinity. The 
post-shock temperature is therefore given by 



where uim is the shock speed in units of 100kms~i, and 
va\o is the Alfven speed in the unshocked gas in units 
of 10kms~i. The thermal and non-thermal spectra are 
matched at an injection energy which is assumed to be 
species dependent: Berezhko & Ellison consider only hy- 
drogen, and specify the injection to thermal momentum 
ratio as A = 4.3. We assume that this momentum rar- 
tio holds for all species, such that the injection energy is 
Einj{j) ~ pf„j/(2mj). By demanding equality of eq.(6) and 
cq(7) at the injection energy, we calculate an injection frac- 
tion of 



77 = 1.85 X 10 •'{m/mj) 



4/3 



(9) 



(8) 



An important consequence of the mechanisms used in the 
present work is that, compared to May ct al. (2000), hydro- 
gen and helium are much more important: in the thermal 
distribution particles have the same energy, not the same 
speed, so impact energy is not weighted to projectiles of 
higher mass. Also, in the non-thermal mechanism, the injec- 
tion energies are higher, and the injection fractions smaller, 
for ions of higher mass, reducing their importance relative 
to light ions. Hence, in contrast with May et al., we cannot 
ignore hydrogen sputtering in the present work. 

3.2 Selective Sputtering of Impurities 

Certain metallic elements, notably Ca, Ti, Co & Ni, are 
observed to be very under-abundant in the gas-phase of the 
ISM (Savage & Sembach 1996). It has long been known that 
gas-phase metallicity correlates with gas velocities (Barlow 
& Silk 1977), suggesting that shocked gas loses some of its 
grain material back to the gas. We test here, the idea that 
such observations may be explained, at least in part, by 
a selective propensity of these elements to resist sputter- 
ing, relative to the bulk grain material, or alternatively, of a 
strong tendency to be stopped, and captured within a grain 
when acting as a projectile. These processes are likely to 
be complementary to the lower energy differential surface 
adsorption proposed by Barlow & Silk (1977) and Barlow 
(1978b). 

The first possibility, a strong resistance to sputtering, 
has been tested by admixing each of the grain materials 
with small percentages of selected elements as impurities be- 
fore carrying out the sputtering calculations. We note that 
some degree of non-stoichiometric sputtering is expected, 
with the more volatile components of the target likely to ex- 
perience a greater degree of sputtering (Tielens et al. 1994). 
Initial calculations were carried out with the SRiM-2000.10 
code, in which the energy parameters. Es, Ed, and Eb are 
not species-specific. These results were negative: that is, the 
fraction of impurity atoms sputtered from the target was not 
substantially different from the initial fractional abundance 
of the same element in the target. In case this negative re- 
sult was due to the uniformity of the energy parameters, we 
also studied selective sputtering from the olivine target us- 
ing the more versatile SRIM-2000.40 which allows separate 
values of the energy parameters to be set for the Mg, Si, 
Fe, and 0-components of the olivine. We note that of these 
consistuents, Mg, Fc and Si are undcrabundant in the gas 
phase relative to oxygen, with Fe being the most extreme 
of the three and Si, the least extreme (Savage & Sembach 
1996). Our results for this case are tabulated in Table 2. 
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Table 2. Total Sputtering Yields for Olivine Constituents 



Element 


% 


% 


ySOOeV 
% 


% 


O 


4.54 X 10"'' 


0.375 


4.010 


3.140 


Mg 


6.24 X 10""' 


0.440 


4.576 


3.761 


Fe 


3.29 X 10""' 


0.224 


2.913 


1.861 


Si 


7.75 X 10"'' 


0.234 


2.810 


1.699 



Notes for Table 2. 

Values for oxygen have been divided by four to allow for its stoi- 
chiometric ratio in the target. 



We show angle-averaged total sputtering yields for each ele- 
ment in the olivine at four different energies. The figures for 
oxygen have been divided by four to allow for its stoichio- 
metric abundance in olivine. The total sputtering yields have 
been calculated by summing over the abundance-weighted 
contributions of the ten commonest Galactic species (see 
Appendix 1). 

The results in Table 2 do show differences between the 
amount of sputtering for each element. Close to the thresh- 
old for sputtering, Si has the highest sputtering yield, which 
reflects its low displacement energy (see Table 1). At higher 
energies we find an approximate binary divide with oxygen 
and magnesium more likely to be sputtered than the iron 
and silicon. This cannot explain the observed sequence of 
gas-phase depletions however: we expect magnesium to be 
intermediate between iron and silicon, which is not the case 
in Table 2, where Mg is more easily sputtered than oxygen 
at all but the lowest energies. 

The second possibility, a strong tendency to be captured 
by a grain, looks far more promising as an explanation for 
gas-phase depletions. The basic argument is that an atom 
from one of the depleted elements, striking a grain as a pro- 
jectile, is likely to be stopped within the grain, and trapped, 
at a depth characteristic of the grain material, and the im- 
pact energy (see Section 3.3). In the event that it does cause 
sputtering, it is likely not to eject an atom of its own kind, 
or one of similar rarity, but an atom of the basic grain mar 
terial: hydrogen, oxygen, carbon, silicon, or perhaps iron. 
We therefore expect a steady loss of such atoms from the 
gas-phase until we reach an equilibrium state, where the 
rate of stopping in grains, for a particular element, becomes 
equal to the rate of return by the dominant grain destruc- 
tion process. Atoms of the heavy elements are more likely 
to be stopped inside the target, with a lower probability 
of backscattering than hydrogen and helium. In Table 3, 
we show the backscattered fraction for projectiles of various 
types, when incident on an olivine target at 0.1 and IkeV. 

The data in Table 3 show that very few ions of the 
elements with large gas-phase depletions get backscattered 
compared with hydrogen, helium and the C,N,0 group of 
elements. Ions from the depleted group (from Ca onwards 
in Table 3) are very likely to be stopped inside the the target, 
and become part of it. The implication is that this is a viable 
mechanism for explaining interstellar depletion, although we 
concede that Ca would not be the most depleted element on 
the the basis of these data. 



Table 3. The backscattered fraction, fg, of ions of various pro- 
jectile ions normally incident on an olivine target at 0.1 and 1 keV 



Element 


/s (100 eV) 


fB (IkeV) 


H 


0.3272 


0.1874 


He 


0.2387 


0.1845 


C 


0.0587 


0.0582 


N 


0.0451 


0.0476 


O 


0.0346 


0.0381 


Ca* 


0.0030 


0.0072 


Ti 


0.0002 


0.0031 


Co 


0.0000 


0.0008 


Ni 


0.0000 


0.0008 


Cr 


0.0001 


0.0019 


Fc 


0.0000 


0.0011 


Mg 


0.0166 


0.0200 


Si 


0.0107 


0.0151 



Notes for Tabic 3. 

* Calcium has the highest gas-phase depletion. 

3.3 Modification of Grain Material 

The 'classical' sputtering yield effectively treats the sput- 
tering process as negative accretion, with a grain radius, a, 
that declines according to the equation 

da _ Yfhnv , , 

dt ~ Ap 

where m is the mean mass of grain atoms, n is the number 
density of projectile atoms in the gas-phase, v is the rela- 
tive velocity of projectile and grain, and p is the density of 
the grain. It is a consequence of the results presented above, 
that this view of sputtering cannot be correct because, given 
that the projectile atoms do not immediately escape the 
material, we find a steady increase in the number of grain 
atoms with time. Wc discuss below whether it is reasonable 
to expect stopped projectiles to escape the grain material 
rapidly compared with other evolutionary timescales. For 
bombarding gas with normal Galactic abundance, most of 
the deposited particles will be light, but we investigate an 
alternative evolution for sputtered grains in supernova rem- 
nants in Section 4, where atoms of the C,N,0 group and 
heavier species have very high abundances relative to typi- 
cal interstellar gas. 

First, we plot the stopping probabilites for the standard 
100 nm grain, with hydrogen projectiles, for the carbon, sili- 
cate and ice materials. The results are shown in Figure 6. We 
can see immediately that apart from a high-energy cut-off, 
dictated by the size of the grain, the fraction of projectile 
atoms stopped within the grain is typically above 0.6. This 
fraction becomes smaller at higher angles of incidence, but 
nevertheless the stopping fraction is, on average, higher than 
the sputtering yield for target materials of all three types. 
The conclusion is that the 'sputtering' process usually re- 
sults in an increase of the number of particles in a grain, 
though not necessarily the grain mass. 

It is vital to know whether embedded projectiles are 
likely to remain embedded within the grain material, or 
whether they can escape on a timescale which is very much 
shorter than the typical duration of the event which intro- 
duced the sputtering. If an embedded projectile can escape 
quickly, then the traditional sputtering yield is effectively 
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equivalent to the true sputtering yield, whilst if the projec- 
tile becomes permanently embedded, the true yield is the 
difference between the traditional yield and the embedding 
fraction; this difference may be negative, leading to grain 
growth. 

Atoms in a solid can diffuse through the lattice, and the 
same mathematical formalism can be used to describe this 
process as is used for gaseous diffusion (Shewmon 1989). The 
diffusion process may be used for both the mobility of im- 
purity atoms, and for 'self-diffusion' of the elements present 
in the lattice. Diffusion through solids is more complex than 
the analogous process in gases because of the 'quantized' na- 
ture of the motion, in which an impurity atom must move 
from one lattice site to another. This extra complexity is ex- 
pressed through a diffusion 'constant' which depends on the 
packing of the lattice, its natural frequency of vibration, the 
temperature of the lattice, and the binding energy of the im- 
purity atom to the lattice. The most important of those do- 
pendencies are the binding energy and temperature, because 
they appear in an exponential Boltzmann factor: it is very 
difficult for a strongly bound atom to diffuse through a very 
cold lattice, as expected. The precise form of the equation 
for the diffusion coefficient depends upon the mobility mech- 
anism: for the vacancy mechanism, where an atom migrates 
by jumping between randomly-placed lattice vacancies, the 
equation is 



Stopping Fraction 



0.8 



0.6 



0.4 



0.2 




2 3 4 

logjd lion Energy (eV)] 



D{T) = aoJ/exp 



R 



exp 



RT 



(11) 



where ao is the lattice spacing, v is & directed oscillation 
frequency (towards a vacancy), R is the gas constant, T is 
the grain temperature, Sx is an entropy and Hx is an en- 
thalpy. The subscript-!/ contributions are from vacancies, 
and the m-subscript contributions refer to lattice activation 
energies. For the interstitial mechanism, where a small atom 
'squeezes' through the lattice by deforming bonds, the equa- 
tion is very similar, but excludes the subscript-j/ terms. A 
diffusion timescale is Tdiff = rf^ / (47r^ (T) ) for spheres of di- 
ameter d (Shewmon 1989). Prom data in Norwick & Burton 
(1975), we can deduce that at 100 K, the measured tempera- 
ture of the processed dust in the SNR of Cas A (Dunne et al. 
2003), all the common elements will be fixed to their lattice 
sites {taudi ff> 10^ yr) except H and He. For example, for 
carbon in iron (assuming the plot in Norwick & Burton can 
be extended down to 100 K at the same gradient) the diffu- 
sion coefficient is of order 10""'* cm^ s"'^, giving a diffusion 
time many orders of magnitude longer than the age of the 
Universe, for even the smallest grains. The fate of the helium 
is not crucial as it has a much lower abundance than H. For 
hydrogen, much depends on whether the H-atoms become 
chemically bound to the lattice. If they do not, for exam- 
ple hydrogen in a- iron (Norwick & Burton 1975), the diffu- 
sion coefficient at 100 K lies between 10~^ and 10~^ cm^ . 
Even assuming the slower value, hydrogen can escape from 
a 10 nm layer on a timescale of 25 /iS. Contrast this with 
the diffusion of H2 through a silicate, where chemical bond- 
ing can take place. Here the diffusion coefficient at 100 K is 
9.8 X 10^^^ cm'^ s^"*^, yielding a timescale of 80000 yr to es- 
cape a 100 nm layer. Hydrogen atoms are more reactive than 
molecules, and so would be expected to remain longer in the 
solid. Chemically bound hydrogen can survive for significant 



Figure 6. The fraction of projectile ions stopped inside the target 
grain as a function of projectile energy for zero angle of incidence. 



times in the solid. Those conclusions are of some importance 
in Section 3.4 and in Section 4. 

Overall, typical deposition rates exceed sputtering 
rates. If the bombarding gas has Galactic abundance, this 
will still lead to net erosion, since sputtering by hydrogen 
and helium, even if these elements remain embedded, re- 
sults in irreplaceable loss of the metallic elements that form 
the grain core. However, in metal-enriched gas, as found in 
supernova remnants, it is possible that embedding of the 
metals can exceed the total rate at which they are sput- 
tered. Hydrogen and helium, if not chemically bound to the 
lattice, can escape from a typical grain on timescales short 
compared to the time needed to erode a grain by sputtering. 



3.4 Processing of an olivine target 

Here, we present results from SRIM computations which ex- 
plicity take into account the modification of a target by the 
combined processes of sputtering and projectile stopping. 
The target material was olivine, as in Section 3.1 and Sec- 
tion 3.2, but modification of the target grain weis studied in 
the following manner: Since SR.im projectiles do not modify 
the target with which they interact, the calculations were 
broken down into a series of steps, each with a fixed tar- 
get composition. Beginning with a pristine target of olivine, 
SRIM calculations were carried out for all ten common projec- 
tile species at various angles and energies. Sputtering yields 
and stopping fractions were recorded, and also the average 
penetration depth for each projectile species, weighted by 
the distribution function developed in Section 3.1. To begin 
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Table 4. Angle-averaged penetration depths for projectile species 
Incident on an olivine target, as a function of energy. At 10 eV, all 
penetration depths were negligible. Columns 2-6 contain penetra- 
tion depths at, respectively, 30, 100, 300, 1000, 3000 and 10000 eV. 
Column 7 contains the energy-averaged penetration depths. 



Element 




dioo 


dsoo 






diok 


d 




nm 


nm 


nm 


nm 


nm 


nm 


nm 


H 


0.7 


1.5 


3.2 


8.4 


20.6 


45.6 


1.746 


He 


0.4 


0.9 


1.9 


5.4 


14.0 


35.2 


1.036 


C 


0.3 


0.4 


0.8 


1.7 


4.0 


11.4 


0.472 


N 


0.3 


0.4 


0.8 


1.6 


3.5 


9.8 


0.472 


O 


0.3 


0.4 


0.7 


1.5 


3.2 


8.7 


0.444 


Ne 


0.3 


0.4 


0.7 


1.4 


2.8 


7.5 


0.444 


Mg 


0.3 


0.4 


0.7 


1.2 


2.6 


6.7 


0.446 


SI 


0.3 


0.4 


0.7 


1.2 


2.4 


5.8 


0.446 


S 


0.3 


0.4 


0.7 


1.2 


2.2 


5.2 


0.446 


Fe 


0.4 


0.5 


0.7 


1.1 


2.0 


4.2 


0.514 



the next step, the target was divided into layers, according 
to the penetration depths, and the composition of each layer 
modified to match a bombardment of 5% by number of the 
original target material. The modification by species was 
carried out in proportion to the sputtering yields and stop- 
ping fractions for that species, and the overall weight of each 
species determined according to Table Al. This use of direct 
abundances, rather than collision rates weighted by thermal 
speeds of the individual species, compensates crudely for 
that fact that the grains in the shocked gas would have a 
large initial drift velocity with respect to the gas, which is 
only slowly dissipated. A fresh set of calculations was then 
performed for the modified target, which now contains pro- 
jectile species. New sputtering yields and stopping fractions 
were computed, and the target modified again to match a 
further 5% by number of the original target. The penetra- 
tion depths, however, were not modified on the second or 
subsequent steps, since stopping distances have very large 
standard deviations, and the number of differentiated target 
layers would rapidly have reached unacceptable proportions. 

3.4- 1 Calculations for the pristine target 

First, we report the penetration depths achieved for each 
projectile type as a function of energy and impact angle. 
Since these values were also used for all subsequent steps, 
we discuss them only once. The depth-data appear in Ta- 
ble 4. The depths were originally measured at energies of 
10, 30, 100, 300, 1000, 3000 and 10000 eV, and angles of 0, 
60, 85°, with the assumption of zero penetration at 90°. In 
the columns of Table 4, the depths have been averaged over 
solid angle, assuming an isotropic distribution of impacts. 
The final column gives the grand-average over all the en- 
ergies used weighted according to the distribution function 
from eq.(6) and cq.(7). 

When integrating over the distribution function, we 
have adopted the following standard parameters: shock 
speed, wioo = 4.0, Alfven speed, vaio = 1.0, and ratio of 
specific heats in the thermal gas, 7g = 5/3. These parame- 
ters are typical of a supernova blast wave, close to its tran- 
sition from the Sedov- Taylor phase to the pressure-driven 
snowplough phase of the remnant. In Table 5, we display 



Table 5. Angle- and energy-averaged sputtering yields for the 
ten common projectile species Incident on a pristine olivine tar- 
get. The final column is the angle- and energy-averaged stopping 
fraction for that projectile. 



Element 






Ysi 


Yo 


fatop 




% 


% 


% 


% 




H 


1.06 


0.73 


0.71 


3.78 


0.435 


He 


5.86 


2.96 


2.90 


19.7 


0.494 


C 


13.2 


7.35 


7.03 


45.5 


0.642 


N 


13.8 


7.68 


7.28 


47.3 


0.661 


O 


14.1 


7.89 


7.44 


48.5 


0.677 


Ne 


14.9 


8.21 


7.71 


50.8 


0.698 


Mg 


15.2 


8.37 


7.79 


51.7 


0.714 


Si 


14.9 


8.10 


7.43 


50.5 


0.735 


S 


14.4 


7.87 


7.06 


49.0 


0.753 


Fe 


13.6 


6.64 


6.35 


45.1 


0.803 



the sputtering results for the pristine target. Each row rep- 
resents a different type of projectile; the first four columns 
are sputtering yields for the four elemental consituents of the 
olivine; the final column is the stopping fraction. The sput- 
tering yields and the stopping fraction are averaged over 
angle and over impact energy, weighted by the distribution 
function. However, the yields are not corrected for either the 
stoichiometric ratio of oxygen in the target, or for the rela- 
tive abundances of the projectile elements. Table 5, together 
with Table 4 and Table Al, provide the means to alter the 
target composition in line with the expected bombardment 
in the post-shock regime. 

3.4.2 Calculations for the modified target 

Following the calculation on the pristine target, the tar- 
get material was divided into layers, dependent on pene- 
tration depth. Each projectile material was assumed to be 
uniformly distributed between the surface and its own pene- 
tration depth, d from Table 4, and to have zero abundance at 
any greater depth in the target. Although, in principle, the 
modified target would have eleven layers (unmodified olivine 
plus ten layers modified by successive mixtures of projectile 
species) the average penetration depth data (final column of 
Table 4) allows a simplification in which the target is repre- 
sented by unmodified olivine plus five modified layers. The 
outermost layer is contaminated by all the projectile species, 
the second by H, He, C, N and Fe, the third by H, He and 
Fe, the fourth by H and He, and the fifth by hydrogen alone. 
Each projectile species was assumed to sputter target atoms 
uniformly from all the zones into which that species could 
penetrate. 

Bombardment of the target was carried out in a se- 
quence of steps, with the first step being bombardment of 
the pristine target. Throughout each step the composition of 
the target, in a given layer, was assumed to be constant. At 
the end of each step, the composition of the target was ad- 
justed to take into account both the loss of sputtered atoms, 
and the embedding of stopped projectiles. At all steps later 
than the first, the projectile species form embedded popula- 
tions, which may themselves be sputtered from the target. 
The very low abundance of all projectiles heavier than He 
made it possible to ignore the sputtering effect of carbon and 
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Table 6. Final abundances in the grain, after 35% bombardment, 
for all penetrated layers, as a fraction of the total number of nuclei 
present in the layer, for the ten most common elements. 



Element 


Layer 1 


Layer 2 


Layer 3 


Layer 4 


Layer 5 


H 


0.1228 


0.1227 


0.1228 


0.1228 


0.1250 


He 


0.0219 


0.0219 


0.0219 


0.0219 




C 


2.41(-4) 


2.24(-4) 








N 


8.0(-5) 


6.8(-5) 








O 


0.4870 


0.4869 


0.4871 


0.4871 


0.4993 


Ne 


1.17(-4) 










Mg 


0.1216 


0.1216 


0.1216 


0.1216 


0.1248 


Si 


0.1231 


0.1232 


0.1233 


0.1234 


0.1254 


S 


1.7(-5) 










Fe 


0.1232 


0.1233 


0.1234 


0.1233 


0.1260 



Layer 1 




Percentage bombardment 



Figure 7. Changes in the abundances of hydrogen and helium 
projectiles and oxygen in a 100 nm grain as a function of total 
bombardment for the outermost layer (see text for definitions). 

heavier elements, even though these heavy nuclei are indi- 
vidually more efficient sputterers then hydrogen and helium 
(see Table 5). Although ignored as sputtering projectiles, the 
embedding of the heavier elements was still considered, and 
populations of these species were also subject to sputtering 
by H and He. For the effects of heavy elements in conditions 
with non-standard abundances (in SNRs) see Section 4. 

The sequence of graphs. Fig. 7 to Fig. 9, shows the 
abundances of all species in the outermost layer of the grain 
as a function of bombardment. The level of bombardment 
is expressed as a percentage, which is 100 x the total num- 
ber of projectile atoms divided by the total number of vul- 
nerable target atoms (summed over all five layers) in the 
pristine target. This degree of bombardment can, of course, 
be converted to time, given a particular astrophysical envi- 
ronment; we consider the effect of supernova shocks below. 
Fig. 7 shows the rising abundances of the sputtering pro- 
jectiles, hydrogen and helium, together with the decrease in 
the fraction of oxygen, the most abundant element in the 
pristine target. Oxygen is the most easily displaced of the 
elements initially present in the olivine. In Fig. 8, we show 
the falling abundances of the metallic components of the 



Layer 1 

0.145 




5 10 15 20 25 30 35 40 
Percentage bombardment 



Figure 8. Changes in the abundances of iron, silicon and mag- 
nesium in a 100 nm grain as a function of total bombardment for 
the outermost layer (see text for definitions). The graphs for iron 
and silicon are so close together that they are virtually indistin- 
guishable. Layer 1 

0.00025 




5 10 15 20 25 30 35 40 

Percentage bombardment 

Figure 9. Changes in the abundances of the 'rare', heavy pro- 
jectiles, C,N,Ne and S in a 100 nm grain as a function of total 
bombardment for the outermost layer (see text for definitions). 



olivine: Fe, Si and Mg. Magnesium is significantly more eas- 
ily sputtered than the other two elements, and we note that 
the sputtering of Si and Fe is so similar that the graphs for 
these two elements are indistinguishable in the figure. Fig. 
9 shows the rising abundances of those projectiles which 
are not considered to cause significant sputtering. The low 
cosmic abundance of these elements means that small con- 
taminant populations develop in the outer regions of the 
grain. 

Overall, the effect of the sputtering process is not so 
much to erode the grain as to process it by significantly 
changing its composition. In all the penetrated zones, we 
find that the total zone mass and the zone density fall below 
the values set for the pristine olivine target, but that the 
particle number in all zones actually rises if we assume that 
H and He projectiles remain embedded in the target. We 
have shown in Section 3.3 that long-term survival of at least 
some hydrogen is likely in a silicate, providing it can bind to 
the lattice. The large number of vacancies produced by the 
sputtering process, leaving a general defHciency of oxygen, 
will aid binding by providing many lattice sites with low 
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energy barriers. The low abundance of He, compared to H, 
makes its loss or retention in the grain of little importance. 



4 SURVIVAL OF GRAINS AGAINST 
SUPERNOVA SHOCKWAVES 

Recent observational studies with SCUBA on the JCMT 
(Dunne et al. 2003) have convincingly demonstrated that 
supernovae (at least those of type II) are net producers of 
dust. We therefore proceed in the spirit of attempting to 
explain this observational fact, rather than attempting to 
make our models agree with previous work, which tends to 
favour net erosion by sputtering. 

The first point that should be made is that a dust grain 
will experience many collisions with gas atoms and ions in a 
typical remnant. In most cases, the grain will receive easily 
enough bombardment to completely erode it if the collisions 
lead to net sputtering over embedding. The mean-free time 
in seconds for collisions between gas atoms/ions and a dust 
grain in the post-shock gas, is given by 



Tcou = 4.087 X 10^ 



^3/8 1/2 
1000 A* 



(12) 



3/4yl/2 2 

where Tiooo is the pre-shock gas temperature in units of 
1000 K, is the relative mass of the bombarding particle 
relative to H, rice is the number density of this bombarding 
species in units of cm~^, vioo is the shock speed in units 
of lOOkms"^, Ts is the post-shock temperature in units of 
10^ K, and dnm is the grain radius in nm. In deriving this 
equation, we have used the definition of the compression 
ratio from Berezhko & Ellison (1999), which is n^/no = 
1.3{vs/cs2)'^^* , whore Vs is the shock speed and Cs2 is the 
sound speed in the pre-shock gcis. For our model grain in 
Section 3.4, being bombarded by hydrogen, the timescale 
for 100% bombardment (being struck by a number of H- 
nuclei equal to the total number of atoms in the vulnerable 
part of the grain) is only 



Tioo = 24.4 



t: 



,3/8 



n t;3/4yl/2 



yr 



(13) 



This figure is largely size independent down to the size where 
the vulnerable zone becomes equal to the entire grain, since 
both the volume of the vulnerable zone and the grain cross- 
section scale as a^. To sputter out all the metallic elements 
from the vulnerable zone requires a few times noo , or about 
a century. For larger grains, including our example, we can 
view grain destruction as requiring the removal of several se- 
quential layers, but oven for our 100 nm grain, total destruc- 
tion would occur on a timescale less than 1000 yr; optimistic 
choices of rice yield ~ 10* yr. The fact that total particle 
number initially rises cannot prevent this destruction: the 
helium, which anyway has a low abundance within the grain, 
would probably diffuse out of the sample in a shorter time 
than Tioo ; hydrogen certainly would if it did not bind chemi- 
cally to the lattice. Although much H might bond, hydrogen 
which did not attach chemically could escape on a very short 
timescale (see Section 3.3). Therefore the true abundance of 
H must always bo loss than in our standard model, and as 
metals are lost, the chances of forming a chemical bond, 
and becoming part of the grain material, are progressively 
reduced. 



In a supernova remnant, abundances are very difi'erent 
to those in our standard model if we are considering material 
which is part of the ejecta, or material which is part of the 
ISM in a young remnant: the early swept material is almost 
certainly stellar wind debris, which contains enhancements 
of the elements He, C, N, and O with respect to the Galactic 
norm. For the ejecta, we have used abundances computed 
by Limongi & Chieffi (2003). Using the Salpeter IMF, the 
mean mass of a supernova progenitor is 18.7 Mq if the upper 
mass limit is set to 100 Mq , but 15.0 Mq if we limit the upper 
progenitor mass to 40 Mq. We have selected the latter figure, 
and abundances, relative to hydrogen, appear in the right- 
hand column of Table Al. To obtain our tabulated values, 
we have summed over all isotopes used by Limongi & Chieffi, 
and averaged over the six models computed for the mass of 
15 Mq. To proceed further, we consider first the fate of the 
supernova ejecta, which are shocked by the reverse shock, 
and then the wind tlcbris and raw ISM, which are shocked 
by the external blast- wave. 



4.1 Dust in the supernova ejecta 

Initially, the ejecta are cold and unshocked, with a contact 
discontinuity separating them from material swept-up by the 
blast- wave (Truelove & McKee 1999). Even when a reverse 
shock is set up by the overpressure of the hot swept mate- 
rial, the cooling times in the ejecta are very short, leading 
to an early radiative phase in the ejecta (Truelove & McKec 
1999). Dust can condense under these conditions within a 
very short time, compared to the evolutionary scale-time of 
the free-expansion phase of the remnant (of order 1000 yr or 
less) . There is observational evidence for dust condensing in 
SN1987A within 530d of the outburst (Dwck 1998). As the 
reverse shock accelerates, higher temperatures are attained 
in the shocked ejecta, but though the bombardment is ef- 
ficient, the ejecta arc metal-rich. Although data in Table 5 
show sputtering exceeding embedding, for metallic elements, 
by a factor of ~ 3, it is quite possible that net embedding 
could predominate. Values of Eb, Ed and Es closer to the 
May et al. (2000) averages would lower sputtering without 
significantly reducing embedding rates. Also, the deposition 
of large amounts of carbon and hydrogen in a thin surface 
layer are likely to lead to the formation of a tough organic 
mantle on grains initially of silicate composition (Greenberg 
& Li 1996) (Li & Greenberg 1997). To compute the sput- 
tering yields from such a mantle, as a function of impact 
energy and time, requires much additional work but such a 
mantle is likely to be at least as resistant as carbon, with a 
yield of about a factor of three lower than for silicates (see 
Figure 2) . The fact that supernovae appear to be net produc- 
ers of dust from observations means that the investigation 
of such mantle formation is worthwhile. Protection of the 
underlying silicate could be effective until typical collision 
energies reach 300-400 eV, allowing the stopping of metallic 
projectiles to exceed their losses in sputtering. Net growth 
of the freshly condensed grains could therefore continue in 
the ejecta until roughly the end of the free-expansion phase. 
Truelove & McKee (1999) quote a characteristic shock speed 
for the reverse shock of. 



vch = 7090^5^ (Mej/Me 



-1/2 



kms 



(14) 
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for supernova energy E^i, in units of 10^^ ergs, with M^j so- 
lar masses of ejecta. The end of the free-expansion phase oc- 
curs when the speed is 0.585 times the characteristic value. 
The mean values of -B51 and M^j for the ISM© model in 
Limongi & CliicfB (2003) arc , respectively, 0.877 and 13.43, 
giving a shock speed of 1060 kms~^. This corresponds, via 
eq.(8), to a temperature of 3.73 x 10® K, and an energy of 
320 eV. The mass of shocked ejecta, at this same stage of evo- 
lution, is 0.781Mej , so most of the ejecta will form dust, may 
accumulate metals, and never bo shocked to a high enough 
energy for net sputtering to dominate. Even the remaining 
22% of the ejecta, which is likely to lose its newly-condensed 
dust to sputtering, will form metal-rich gas which, when it 
eventually cools, is likely to add to the dust already present 
in the cooler shocked ejecta. We note that the cooling time 
is given by 

j.3/2 

Tcooi = 6.3 X 10"^- yr (15) 

where rice is the number density of the cooling gas in cm^'' 
and Cm is a metallicity correction factor (Ciofii, McKee & 
Bertschinger 1988). Therefore, the destruction of the dust 
enhances C,m, and enhances cooling towards temperatures 
where the bombardment process results in net grain growth. 
Overall then, much of the metal content of the ejecta can 
survive the supernova remnant bound in dust. We note that 
the survival of initially condensed ejecta dust was predicted 
by Barlow & Silk (Barlow & Silk 1977), but this would pro- 
vide insufficient mass to make supernovae net creators of 
dust. 

4.2 Dust in the ISM and wind debris 

The fate of grains in the wind debris and raw interstellar 
gas is more complicated. We expect the wind material to 
have enhanced abundances (over the mean ISM values) of 
lower-mass rmclcosyrithesis products, but not to have high 
abundances of the iron-group elements. We also expect suc- 
cessive shells of wind debris to have come from a sequence 
of evolutionary stages with different chemical composition. 
Nevertheless, we expect much of this material to resemble 
red supergiant winds, and to be rich in dust. For most rea- 
sonable parameter values, the entire wind-zone will be swept 
within the free-expansion phase of the remnant, and will 
therefore be shocked and heated to very high temperatures. 
For the particular case of the mean Limongi & Chiefi (2003) 
15 M© model, the blast-wave speed at the end of the free- 
expansion phase is still 1665 km s~^ which is destructive, 
even for metal-rich gas. However, as in the inner regions of 
the ejecta (see above) the dust destruction leads to metal- 
rich gas with cooling times that are much more rapid than 
those for typical ISM gas. With the help of the Berezhko 
& Ellison (1999) formula, n^/no = 1.3(vs/cs2)^^'*, we can 
calculate a compression factor of 113 for the end of the 
free-expansion phase, and a corresponding temperature of 
6.1 X 10® K from eq,(8). With typical coolants at levels of 
100 times their typical ISM abundances, eq.(15) reduces to 
84/ncc yr which is of similar order to the grain destruction 
time. It is therefore likely that the shell of wind debris be- 
comes rapidly radiative, and that dust either survives or 
rapidly recondenses, joining the mass of dust formed in the 
ejecta. This mass of dust is likely to be similar to, but not 



greater than, the mass of dust in the ejecta, on the grounds 
that the pre-supernova star was unlikely to have lost more 
than half its ZAMS mass as wind debris. 

Once the blast-wave heis cleared the wind material, usu- 
ally before the end of the free-expansion phase of the rem- 
nant, it begins to sweep-up unenriched ISM gas, which is 
likely to contain dust at typical ISM proportions of ~1% 
by mass. This gas is still heated to high temperatures, and 
now there is no chance of net embedding for elements heav- 
ier than H and He (sec Section 3.4). Gas bombardment of 
grains in this regime can only be destructive, and nor will 
the release of metals into the gas reduce the cooling time 
very much. This state of affairs will last until the end of 
the Sedov- Taylor phase, that is until the shock becomes ra- 
diative, and cooling times, even for unenriched gas, become 
short compared to the evolutionary timescale of the rem- 
nant. Some dust destruction may continue into the early 
PDS (pressure-driven snowplough) phase which follows the 
Sedov- Taylor era. 

A rough estimate of the amount of dust destroyed can 
be obtained from the mass swept by the blast-wave dur- 
ing the Sedov-Taylor (ST) phase of the SNR. The PDS 
phase begins after ~ 13000n^^*''^("m^^^'* yr (Cioffi, McKee 
& Bertschinger 1988). From the end of the free-expansion 
phase up to this time, the remnant sweeps up a mass of 

M = 17. 5 fl{Mej/MQ) Mq (16) 

which is notably independent of the ISM density. The quan- 
tity p, is the mean relative particle mass in the unshocked 
ISM. For our standard model at 15 M©, we find a total swept 
mass of 306 Mg , so for a typical dust mass fraction about 
3 M© of pre-supernova dust is destroyed. 



5 DISCUSSION 

The results of this work are broadly as follows: Firstly, sput- 
tering in the classical sense is almost invariably an inefficient 
process: In most situations, the sputtering yield for a hy- 
drogen projectile, and the chosen targets, is far lower than 
the 100% required to ensure erosion of the grain. Slow ero- 
sion is still possible at sputtering yields below 100%, but 
this requires that most of the incident projectiles do not get 
stopped in the target and hence become part of the grain. 
In fact, we find that the stopping fraction is typically large, 
and ususally, the number of atoms in the grain will increase 
due to particle collisions, as a result of the stopping of pro- 
jectile atoms within the grain. In many cases, the grain mass 
will actually increase. For particle collisions to be effective 
at destroying grains, we need to appeal to the secondary 
effects of lightening the average atomic mass in the grain, 
the disruption of its crystal structure and diffusion of light 
elements from the grain. 

The analysis of survival of dust in a supernova remnant 
produces an ambiguous result, but we cannot exclude net 
production on the basis of this very simple estimate. The 
Sedov-Taylor phase destroys ~ 3M© of ISM dust. On the 
other side of the balance sheet, the mass of elements heavier 
than He in the ejecta is ~ 1.2 M©, so the ejecta and wind 
debris could between them contribute about 2M© of dust. 
The masses of dust produced and destroyed are very similar, 
and it is necessary to complete a very detailed analysis to 
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find out which predominates, though it appears to be a re- 
quirement from models of galactic dust evolution (Edmunds 
2001) that supernovafi must be net producers of dust. 

The dust produced by condensation followed by bom- 
bardment would have an interesting structure: the core 
would be a condensate, possibly crystalline initially, but 
bombardment would then insert heavy elements at random 
down to their typical penetration depths in the grain. Thus 
even if the original condensate were a silicate, the outer lay- 
ers would be rich in heavy elements that did not form part of 
the mineral: carbon and nitrogen for example. This forma- 
tion and processing scheme could make grains very similar 
to the silicate core with 'organic' mantle type that form the 
large grain population in Greenbcrg's unified model (Green- 
berg & Li 1996) (Li & Greenberg 1997). The dust is also 
formed in the correct region, as observationally determined: 
between the blast- wave and reverse shocks (Dwek 1998). 
Changes in dust grain structure due to bombardment could 
also affect the optical properties of the grains, and this may 
be required to explain the SCUBA observations of supernova 
dust emission compared with more general ISM dust. 
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APPENDIX A: ABUNDANCES 

In our standard model, projectile ions were assumed to have 
Galactic abundances, and the ten most common elements 
were used, with abundances, relative to hydrogen, listed 
in the second column of Table Al. All other elements, at 
Galactic abundance, have individual abundances lower than 
3.58 X IQ-" relative to hydrogen. 

The right-hand column of Table Al shows abundances 
in supernova ejecta computed by Limongi & Chieffi (2003) 
for a 15 Mq progenitor. These figures are the sums over all 
isotopes shown, and averaged over the six models run for 
this mass. 



